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We report on some technical aspects of our calculation of corrections to R{s) and the semi-leptonic r decay 
width [Il2l3j . We discuss the inner structure of the result as well as the issue of its correctness. We demonstrate 
recently appeared independent evidence positively testing one of two components of our full result. 



1. Introduction 

Three important physical observables, namely, 
the ratio R(s) ~ "^^f / -*hadrons) ^j^^ hadronic 
decay rate of the Z-boson and the semileptonic 
branching ratio of the r-lepton are expressed 
through the vector and axial-vector current cor- 
relators (see, e. g. reviews |4I5) ). Perturbative 
QCD provides reliable predictions for these cor- 
relators in the continuum, i.e. sufficiently above 
the respective quark threshold and the respective 
resonance region. 

The 0{a^) result for the massless vector corre- 
latoiH has been known since many years |6l7j . 

Recently the calculation of the next, order 
contribution to the vector correlator has been per- 
formed [l(2j3j . The aim of the present work is 
to discuss some technical aspects of our calcula- 
tions as well as to provide some new arguments 
in favour of their correctness. 

Due to lack of space no phenomenological im- 
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work we will discuss non-singlet contributions only. 



plications of [112] are discussed and the interested 
reader is referred to Refs. [1181911011111^ . 

2. Generalities 

Consider the two-point correlator of vector 
quark currents and the corresponding vacuum po- 
larization function (j^ — Qjfj,Q; Q is a quark field 
with mass m, all other Uf ^1 quarks are assumed 
to be massless) 

= {-g^.uq^ + q^.qu)^{q^)■ (1) 

The physical observable R{s) is related to I\-{q^) 
by 

R{s) = 12Tx^n{q^ + (2) 

For future reference it is convenient to decom- 
pose R{s) into the massless contribution and the 
one quadratic in the quark mass as follows (a^ = 
as{^l'^)/^:)■. 

R{s) = 3|ro^-H^r,^} + ... 




1 



2 



The corresponding representation for the polar- 
ization function reads (Q^ = —q^) 

2 

n = no{L, as) + ^n2(L, as) + Oil/Q^). (3) 
Note that both function on the rhs of ([3]) depend 

2 

on only as and L = In^ and could be conve- 
niently decomposed as follows (n = 0, 2) 

n„ = ^nra:, nr= ^ n^,L^. (4) 

i>0 0<j<i+l 

The terms in (j4]) without L-dependence do not 
contribute to R{s). 

For the calculation of r^'"^ the divergent parts 
of five-loop and the finite parts of the four-loop 
diagrams are needed [.13^. The organization of the 
calculation is best based on using of the evolution 
equation for 11 (see, e. g. [14 1) 



±Uo = j''''ias)-Uas)as-^]llo, (5) 



where = J2i>o'yY^ '^l (subtractive) 

anomalous dimension of the correlator ([1]) and 

/3 = — J2i>o Pi c'i*^^'' is the QCD /3-function. 

To evaluate the L-dependent pieces of the 
polarization function IIq . . . Iljj in terms of 
7^^ . . . 7]"^ and no . . . n° the evolution eq. © 
can be solved perturbatively: 



Tll=^rL + n%^ U\ = ^y^L + Tl%, (6) 

= /3o ^ + i {iV + /3o n%) -t- (7) 



+ ^ 71^^ + 2/30 72^^ + 2/30^0%) 

+ L + n?^o + 2 /3o n§^o) + n§,o, 



(8) 



+1 (3a Pi n% + 3 /?2 n§^o) + L (71-^ + [i^ n% 
-i-2/3ino,o + 3/5on%) + n%. 



(9) 



The evolution equation for 112 describing the 
m^-corrections looks similar to ([5|), namely |3|15j . 

^n2 = -(^2 7™(a.)+/3(a.)a.^^n2, (10) 

where 7m = ^ S j>o '^™^^^'^s is the quark mass 
anomalous dimension. 



3. Results 

We refer the reader to [1] for a discussion of 
various theoretical tools used to compute IIq to 
order and 7^^ to a^. We only want to mention 
here the indispensable role of the parallel version 
[Hj of FORM [17] and the availability of large 
computing resources. The results for 7^^ and 11 
are given in the next two subsections. 



3.1. Five loop anomalous dimension 
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3.2. 0{al) polarization operator 
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12 ' 24 9 
The results for the analytical calculation of 
n§ have been reported in [3] while the 
four-loop quark anomalous dimension is known 
from [18119] and three-loop QCD /3-function from 
[20121) . Numerically we find for the polarization 
operator 



0.0422172 + 0.0253303 L, 
-0.00569664 + 0.0253303 L, 
0.0457538- 0.0080559 n/, 
-l-(0.0502986- 0.00292047 n/)L, 
+ (0.0348292 - 0.00211086 n/) L^, 



^^3 



(18) 
(19) 
(20) 



(21) 



+0.23570 - 0.033603 Uf + 0.0007394 nj 
+ (0.462093 - 0.10679 n/ + 0.0021836 n^) L 
+ (0.219061 - 0.02644?!/ + 0.00048674 n^) 
+ (0.063854 - 0.0077398 n/ + 0.0002345 nj) 
and 



= 

n? = 
nl = 



Hi = 



-0.151982, 

-0.405285 - 0.303964L, 
-4.27066 + 0.200532 n/, 
+ (-3.20428 + 0.109765n/) L, 
+ (-0.721913+ 0.0253303 n/) 



(22) 
(23) 
(24) 



(25) 



-53.0381 + 5.21239 Uf - 0.0740141 nj 

+ (-43.9568 + 4.05526 n/ - 0.0586349 nj) L 

+ (-14.2641 + 1.1082 n/ - 0.0182941 n^) 

+ (-1.80478 + 0.143538 n/ - 0.00281448 nj) L^. 

Note that at eqs. pim3ll5l[TC)) and P^20I221|M| 
as well as L-dependent pieces of of eqs. (|21l25p 
are, in fact, known since long [6I7I22I15] . 
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Table 1 

Possible irrational structures which are allowed 
to appear in £-loop massless propagators. 
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Table 2 

Possible irrational structures which are allowed 
to appear in i-looj) anomalous dimensions and f3- 
functions. 



4. Final results for R{s) 

Our final results for r^'* and r^'* are easily 
obtained from results listed in the previous two 
subsections. Explicit expressions can be found in 

It is of interest to discuss the structure of tran- 
cendentalities appearing in eqs. (|imi7p . On gen- 
eral grounds one could expect that the variety 
of (^-constants entering into MS-renormalized (eu- 
clidian) massless propagatorsQ should depend on 
the loop order according to Table [TJ Table [2] pro- 
vides the same information about possible irra- 
tional numbers which could show up in anomalous 
dimensions. Table [2] comes directly from Table [T] 
by noting that any £ + 1-loop anomalous dimen- 
sion can be obtained from properly chosen £-loop 
massless propagators [23]. 

An examination of eqs. (|imi7p immediately re- 
veals that the real pattern of trancendentalities is 
significantly more limited than what is allowed by 
Tables [T] and [21 Indeed, the four-loop anomalous 
dimension 'y^^ contains no ^4 and no (5 while 
the three-loop polarization operator contains (5 
but does not comprises ^4. The fact of absence 



^It is understood that 0(e(5~*) terms in an £-loop massless 
propagator could contribute only to six-loop anomalous 
dimension and, thus, are not constrained by Table [T] 
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of (4 in 0{al) contribution to the Adler function 
is well-known and well-understood |6l24j . Why 
73^^ is free from C5 is still unclear (at least for 
us). 

Let us move up one loop. The situation is get- 
ting even more intriguing: the five-loop anoma- 
lous dimension does contain (^4 but still does 
not include CliCe and ^7. The four- loop polar- 
ization operator contains ^4 but is free from ^g. 
Even more, after we combine and II'^ to pro- 
duce the Adler function, the resulting coefficient 
in front of C4 happens to be zero in a non-trivial 
way! Indeed, the contribution proportional to ^4 
from n§ Q reads 

(55 55 ^ \ _ ( 605 I 715 „ 55 ^2 A 

and is exactly opposite in sign to the correspond- 
ing piece in (fT4|)! 

Unfortunately, we are not aware about exis- 
tence of any ratio behind these remarkable ob- 
servations. 

5. How reliable are our results? 

The history of multiloop calculations teaches 
us to be cautious. For instance, approximately 
twenty years ago a severely wrong result for the 
0(q;^) coefficient in R{s) was pubhshed [25] and 
corrected only three years later [6|7j . 

Now one of these authors (rightfully!) rises an 
important issue of the correctness of the results 
[ll2j and emphasizes the necessity of performing 
their independent test |26|27j . 

We completely agree with this argumentation. 
Unfortunately, at the moment, we are not aware 
of any independent team which is going or, at 
least, able to check our results in full. 

However, as described below, the results of a 
recent calculation [28] allow at least for a (partial) 

test of [nas]. 

5.1. A test of the polarization operator 

In Ref. [28j a large amount of information 
about the massive four-loop polarization function 
was collected (its threshold behavior [29|30|31] . 
as well as low-energy moments [32133134] and 
high-energy asymptotic *15l35j) in order to re- 
store the whole function within the Padc ap- 



proach [36|37|38j by properly extending the treat- 
ment elaborated more than a decade ago for the 
massive three-loop polarization Junction 39J . 

Within this method it is customary to deal with 
a "physically" normalized polarization operator 
n defined such that 

n(M, Q, as) = n(M, Q, a,) - n(Af, Q = 0, a,), 

where n(Af, Q,as) is defined by eq.([3|) with the 
use of the pole quark mass M mass instead of 
MS renormalized quark mass m (see [40141142] '). 
Using the results for n(Af, Q = 0, as) as listed 
in [32|33| we arrive at the following asymptotic 
behavior of the (four-loop part of) n(M, Q, Os) 
at Q — * oo. 

fto = (_Q 0638534543 (26) 
-h0.007739812639345n/ - 0.0002345397769 nj) 

(0.219061347 - 0.0264409511 n/ 
-hO.000486744424 n^) 

+L (-0.4620927910 0.1067886396 n/ 

0.0021836455422 n)) + H^J^^ , 

= -11.4121461108 (27) 
H-1.4413529302n/ - 0.032814657849 n^, 

-I = —(-0.4511958959073 (28) 

Q-^ z 

-1-0.03588458587 Uf - 0.00070361933 n}) 

+— (3.084755098861 

z 

-0.2601632475816 nf + 0.004573525651 nj) 

+ - (-6.6515904245 

z 

-H0.78236916962n/ - 0.01465873606 n^) 



H\^> = -8.16060818463927 (29) 
-hi. 0812664904869 n/ - 0.031095026978 n^. 

To be in agreement with the notations of [55] 
we have used in (|26l28p L = \ii{Q'^ / AP) , z = 
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71 f - 4 


Uf ^ 5 




Fade 


-6.122 ±0.054 


-4.989 ±0.053 


o-(3) 
^0 


exact 


-6.17176892 


-5.02574791 




Fade 


-3.885 ±0.417 


-3.180 ±0.405 




exact 


-4.33306265 


-3.53165141 



Table 3 

Comparison of the Fade method predictions for 



Hq^'' and H\^' with exact results. 



r(3) 



~Tm^ and set the renormalization scale fi — M. 
In addition, we have separated in eqs. (|26l28p in- 
dependent pieces (known since long) from the new 
L- independent ones. 

The authors of [55] did not have at their dis- 

(3) (3) 

posal the L-independent terms Hq and 
and, thus, did not use them. Instead, they were 
able to reconstruct these terms for two partic- 
ular values of n/ = 4,5 from their final Fade 
approximants. Table [3] compares their (approxi- 
mate) results with our exact ones. We observe the 
full agreement (within accuracy of the Fade ap- 
proach) between our (exact) results and approxi- 
mate ones obtained in [28] . 

5.2. Discussion 

The full result for i?(s) is composed from two 
parts: the four-loop polarization operator (in- 
cluding its constant, that is ln((5^) independent 
terms) and the five-loop anomalous dimension 
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At the level of separate Feynman diagrams 



the evaluation of both parts is reduced to the di- 
rect calculation of four-loop massless propagators 
[l4] . On the other hand, the input data used in 
[28] came from three different sources: 

(i) the threshold behavior of the polarization op- 
erator; 

(ii) the ln((5^) dependent part of the high energy 
limit of the four-loop polarization operator which 
technically was obtained by a calculation of three- 
loop massless propagators only; 

(iii) the first two physical moments of the four- 
loop polarization operator which technically were 



obtained by a calculation of massive tadpoles only. 

Thus, we consider the full agreement demon- 
strated by Table[3]as a non-trivial and completely 
independent confirmation of the correctness of 
the results of [1.2.3) . 
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